Dynamics of Black Hole Pairs II: 
Spherical Orbits and the Homoclinic Limit of Zoom-Whirliness 



OO 
O 

o 

(N 

> 
O 

(N 



Rebecca Grossman** and Janna Levin*'' 
** Physics Department, Columbia University, New York, NY 10027 
'Department of Physics and Astronomy, Barnard College of Columbia University, 3009 Broadway, New York, NY 10027 
Institute for Strings, Cosmology and Astroparticle Physics, Columbia University, New York, NY 10027 

becky (Sphys.columbia.edu and 
janna@astro. Columbia, edu 

Spinning black hole pairs exhibit a range of complicated dynamical behaviors. An interest in 
eccentric and zoom- whirl orbits has ironically inspired the focus of this paper: the constant radius 
orbits. When black hole spins are misaligned, the constant radius orbits are not circles but rather 
lie on the surface of a sphere and have acquired the name "spherical orbits". The spherical orbits 
are significant as they energetically frame the distribution of all orbits. In addition, each unstable 
spherical orbit is asymptotically approached by an orbit that whirls an infinite number of times, 
known as a homoclinic orbit. A homoclinic trajectory is an infinite whirl limit of the zoom- whirl 
spectrum and has a further significance as the separatrix between inspiral and plunge for eccentric 
orbits. We work in the context of two spinning black holes of comparable mass as described in the 
3PN Hamiltonian with spin-orbit coupling included. As such, the results could provide a testing 
ground of the accuracy of the PN expansion. Further, the spherical orbits could provide useful 
initial data for numerical relativity. Finally, we comment that the spinning black hole pairs should 
give way to chaos around the homoclinic orbit when spin-spin coupling is incorporated. 
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A complete knowledge of the dynamics of black hole 
pairs is essential for future gravitational wave experi- 
ments. Yet the importance of dynamics has not always 
been appreciated. Although stellar mass black hole pairs 
are significant candidates for a first direct detection with 
LIGO, their detectable gravitational radiation would be 
emitted from nearly circular orbits - at least that was the 
refrain. This preferential focus on quasi-circular inspiral 
was motivated by considerations of long-lived binaries 
that begin with a modest eccentricity that is gradually 
shed as angular momentum is lost to gravitational waves. 
A fair assessment of known astrophysics, the claim dis- 
couraged analyses of orbital dynamics in favor of the sim- 
pler analysis of circular orbits. 

However, black hole binaries formed by tidal capture 
in dense star clusters, such as globular clusters, do not 
conform to this story [l| . As one black hole scatters with 
another black hole in a dense region, a burst of radiation 
is emitted on close encounter. Some subset of these en- 
counters will leave the pair bound in a highly eccentric 
orbit that merges too quickly to circularize. Estimates 
conclude that as many as 30% of multi-black hole sys- 
tems will retain eccentricities > 0.1 as their waves sweep 
into the LIGO bandwidth 0. 

Most recently, a new source of eccentric mergers was 
predicted to have a substantial detection rate [3]. Black 
hole/black hole scattering in galactic nuclei would sim- 
ilarly lead to tidal captures and highly eccentric, short- 
lived black hole binaries with 90% entering the LIGO 
bandwidth with eccentricities > 0.9 Q- These compet- 
itive sources for a first detection by Advanced LIGO 
further motivate our study of the complete dynamics of 
binary black holes 0, H, Q. 

In paper I in this series 0, we presented the spec- 
trum of orbits in the strong-field regime when only one 



body spins. ^ There we found zoom- whirl behavior - dur- 
ing which an orbit zooms out in large leaves followed by 
nearly circular inner whirls. Significantly, this extreme 
form of perihelion precession is prevalent in comparable 
mass systems 0, just as it is in extreme-mass-ratio in- 
spirals 0]. Zoom- whirl behavior is not restricted to ex- 
treme eccentricities, but can be executed by orbits of all 
eccentricities in the strong-field. We should therefore be 
prepared to detect evidence of such black hole dynamics 
in gravitational waves. 

In this companion to paper I, we work again in the 
conservative Hamiltonian 3PN approximation plus spin- 
orbit coupling, but move beyond paper I to consider two 
spinning black holes in a binary. We focus on special sets 
of orbits, namely the spherical orbits. That might seem 
ironic since we have just argued that gravitational wave 
science will probe the full range of dynamical possibili- 
ties, there are several good reasons to devote some time 
to constant radius orbits. 

1. If even one black hole spins, the constant radius or- 
bits are no longer circles (unless spins are exactly 
aligned or anti-aligned with the orbital angular mo- 
mentum) [T^]. As a result of spin precession, 
they fill a band on the surface of a sphere and have 
thereby been coined spherical orbits (see Fig. [1]). 
Long-lived binaries that have shed enough angu- 
lar momentum to lose their eccentricity but not 
their spin will exhibit quasi-spherical inspiral and 
not quasi-circular inspiral. To detect waves from 
realistic binaries, we will need to understand the 



^ In this paper we will use the word "orbit" to mean bound, non- 
plunging trajectories. 
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FIG. 1: A spherical orbit for mass ratio milmx — 1/4 and 
spin amplitudes of 3/4. Both spins are initially displaced from 
the orbital angular momentum by 7r/4. Notice the orbit is not 
closed. Upper Panel: The three-dimensional orbit fills out a 
strip on a sphere. If we waited long enough, the band would 
be solidly painted, a reflection of the aperiodicity of the orbit. 
Lower Panel: The path as caught by the orbital plane reveals 
the constant radius. 

orbital parameters and precessions of these spheri- 
cal orbits. 

2. Black hole pairs that enter the LIGO bandwidth 
with their eccentricity intact will evolve through a 
sequence of zoom-whirl orbits rather than nearly 
spherical ones. Still, the spherical sets are special 
since they mark the minimum and maximum en- 
ergy in the strong-field spectrum of bound orbis 
for a given angular momentum.^ The orbital de- 
mographic is therefore entirely determined by the 
spherical orbits . 

3. The energetically-bound unstable spherical orbits 
mark the divide between inspiral and plunge 
More specifically, black hole spacetimes harbor ho- 
moclinic orbits - orbits that approach the unstable 
spherical orbits in the infinite future or in the infi- 
nite past [1, [13, Elj [3 • ^ Homoclinic orbit, often 



refered to as a separatrix, is a classic feature of 
a non-linear dynamical system and deserve atten- 
tion since they mark the transition from inspiral to 
plunge for all pairs. As we mention in the close 
of this paper, when spin-spin coupling is incorpo- 
rated in the PN Hamiltonian, the homoclinic set 
can become the locus of a transition from regular 
to chaotic behavior. 

Among the infinite list of spherical orbits, two are 
valuable enough to deserve names: the innermost stable 
spherical orbit (isso) and the innermost bound spherical 
orbit (ibso). The acronyms are drawn in analogy with 
the equatorial isco (innermost stable circular orbit) and 
ibco (innermost bound circular orbits). The isso is the 
lowest energy spherical orbital and the ibso is the highest 
energy, bound spherical orbit. 

The isco, well-known as the site of the transition from 
inspiral to plunge for quasi-circular inspiral, is actually 
the zero eccentricity homoclinic orbit [l^ . All other 
orbits, besides the quasi-circular one, will transition to 
plunge through another member of the homoclinic fam- 
ily, hence the importance of the homoclinic set to grav- 
itational wave science [l^, To our knowledge the 
homoclinic orbits have not yet been identified in the PN 
Hamiltonian expansion before this paper, although an 
earlier paper found homoclinic orbits and zoom-whirl be- 
havior in a hybrid PN expansion j^]. Excitingly enough, 
homoclinic orbits have been observed in fully relativistic, 
numerical treatments as well [T5| . 

Taken together this special set - composed of spheri- 
cal and homoclinic orbits - demarcates dynamical regions 
and we spend time on their attributes in this paper. Due 
to the lack of confidence in the PN expansion at close sep- 
arations, we do not advocate that these results be taken 
as quantitatively accurate descriptions of binary black 
hole dynamics, but rather as qualitatively descriptive."^ 
Inecd, we point out peculiar artifacts of the 3PN system 
as we go along and the results of this paper could provide 
a new terrain on which to test the PN expansion against, 
for instance, numerical relativity. 

Our approach has some overlap with, but is not re- 
dundant with, the Refs. 0, |T3| and allows us to find 
homoclinic orbits and stability exponents for use in the 
periodic orbit taxonomy of paper I Q ■ We also simplify 
the initial conditions for spherical orbits in the absence 
of radiation reaction. For quasi-spherical orbits with ra- 
diation reaction included see [l0| . 

The outline of the paper is as follows: In ^ we write 
out the equations of motion for two spinning bodies in an 
orbital basis, relying on the results of appendix [XI In f|TT] 
we determine the orbital parameters of spherical orbits. 
In jjllll we find the homoclinic orbits and emphasize their 



^ As detailed in the paper, we define the strong-field by the ap- 
pearance of bound, unstable spherical orbits. 



^ The weakness of the PN approximation famously plagues other 
attempts to pinpoint the transition from inspiral to plunge 
through the isso llC| . 
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connection to dynamical instability. In the conclusion, 
mV[ we discuss the destruction of the spherical orbits 
and the transition to chaos when spin-spin coupling is 
included. 



I. 3PN HAMILTONIAN + SO COUPLING 

We will work with a condensed and revealing set of 
equations of motion in a non-orthogonal orbital coor- 
dinate system as derived in For reference in this 
companion to that paper, we write out the usual 3PN 
Hamilton plus spin-orbit coupling for two spinning black 
holes. 

In a Hamiltonian formulation, the equations of motion 
are derived from 



dH 
dp 



dH 
dr 



(1) 



As is standard convention, we work in dimensionless co- 
ordinates: the dimensionless coordinate vector, r, is mea- 
sured in units of total mass, AI = mi + m2, for a pair 
with black hole masses mi and 1712- The canonical mo- 
mentum, p, is measured in units of the reduced mass, 
fi = mim2/M . The dimensionless combination = fJ-fM 
will prove useful. We write vector quantities in bold. 
The coordinate r is to be understood as the magnitude 
r = y/r ■ r. Unit vectors such as n = r/r will addition- 
ally carry a hat. Finally, we have used the dimensionless 
reduced Hamiltonian H = H/fi in Eqs. ([T]), where H is 
the physical Hamiltonian, to 3PN order plus spin-orbit 
terms 0, [H, S [Si ■ H can be expanded as 



H — Hn + HiPN + H2PN + H-iPjy + Hso 



(2) 



where 



HiPN = ^ (3^ - 1) (p')' - ^ [(3 + V) P' + '?(ft • P)'] \ + ^ 
H2PN^^{1-57j + 5tj^) {P^f + l 
-2r;2(n.p)V-3r7^(n.pr]^ 



(5 - 2077 - 37^2) (p2) 



27- 
2 



+ ^ [(5 + 87y) p2 + 37y(n • p)^] 1 - i (1 + 37y) 1 



H3PN = (-5 + 357; - 7O772 + 35773) (p2)^ -I- 1 ^ 



+ (2 - 3ii)if{h ■ p)'{p'f + 3(1 - 7j)7j'{h ■ p)^p' - 5773(6 • pf] 



-7 + 4277-53772 - 577^) (p2)' 
1 

r 



-^{-27 + 13677 + 109772)(p2)2 + -l(i7 + 3077)77(11 • pfp^ + ^(5 -f- 437;)7;(fi • p)" 



16 



16^ 



12 



{ l^ ^^^^^ ^ ^^^^ " ^^"^^^ ^' ~ 55271"^] P^ - ^ (340 + 37r2 + lUij) rj{h ■ pf 



1 



-[12 +(872 - 63^2)^1 



H 



96 

L • Seff 



SO 



r 



(3) 



(4) 



For two spinning black holes Sgff is^ 
Seff = <^iSi + ^282 



(5) 



* The definitions for Seff can vary in the literature up to an overall 
constant although the reduced Hgo must be the same for all 
prescriptions. 



where the dimensionless reduced spins are defined as 

Si =ai(m2/^Af) , S2 =a2(7772/^Af) . (6) 
and 



^1= 2 



3m2 
2mi 



77 



3mi 
2m2 



(7) 



The dimensionless spin amplitudes are confined to the 
range < ai.2 < 1- The reduced orbital angular momen- 
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turn L = r X p and the spins precess according to 

Seff X L 



L 

Si 
S2 



= s- 



r3 

L X Si 
L X S2 



The spin precessions can be grouped together, 

. . L X Seff 

+ f 2 — 5 • 



(8) 



(9) 



Notice that the precession of the sum of the spins is equal 
and opposite to the precession of the orbital angular mo- 
mentum. So that J = L-|-Si-|-S2is conserved. The 
magnitudes L, and S2, the inner product Seff • L, and 
the energy (the Hamiltonian) are also constant for a given 
orbit. 

In general, neither J -L nor the magnitude \Son\ is con- 
stant. However, there arc notable exceptions [24I. Both 
J • L and the magnitude l^eff | are constant (1) if one of 
the black holes is spinless as was the case in [7|, (2) if the 
binaries have exactly equal mass [2^ . or (3) if both spins 
are aligned or anti-aligned with the angular momentum. 
Case (1) is worked out thoroughly in paper I 0]. To see 
that the claim is true in the equal mass case (2), notice 
that Seff = 5i(Si -I- S2) and therefore J = L + Seff/(5i. 
Consequently, J • L = + Seff ■ L/5i is conserved since 
both terms on the right hand side are conserved. Fur- 
thermore, Soft = <^i(Si -I- S2) and it follows from Eq. ([9]) 
that the change in Seff is always perpendicular to Seff 
so its magnitude remains constant. In case (3), when the 
spins are aligned or anti-aligned with the orbital angular 
momentum, motion is confined to a plane and there is no 
precession. Therefore Seff and L are constants and the 
rest follows. 

The results of this paper will apply to a general Seff 
unless explicitly stated otherwise. 



A. Equations of Motion in the Orbital Plane 

In a non-orthogonal orbital basis, the equations of mo- 
tion assume a simple form that allows us to analyze the 
dynamics of the black hole pairs. The plane perpen- 
dicular to the precessing orbital angular momentum is 
spanned by the vectors (ii, €>) where n — r/r and 



^' = L X fi 



(10) 



Notice n • ^ = so these basis vectors are orthogonal. 
The entire orbital plane then precesses around the con- 
stant total angular momentum J = L + Si + S2 in the 
direction ^ defined through 



* = J X 



(J 


X L) 


J 


X L 



(11) 




JxUl-0) 




Orbital plane 



FIG. 2: Upper: The orbital plane precesses around the J = k 
axis through the angle 'P. Lower: The orbital plane can be 
spanned by the vectors (n, 4>). 



The construction, familiar from classical celestial me- 
chanics [H, m, [2i, i^, is illustrated in Fig. H Inci- 
dentally, this basis is explicitly constructed for J x L 7^ 
0. When the spin and orbital angular momentum are 
aligned, anti-aligned, or spin is zero then J x L = 0, mo- 
tion is confined to a plane, and we should use the usual 
equatorial planar basis. 

Notice that ^ is not orthogonal to the orbital plane 
and therefore our orbital basis (fi, ^,4') is not orthogo- 
nal. 

From Eq. ([1]), we can find equations of motion 
in coordinates (r, <&, 5*) and their canonical momenta 
{Pr , P<t, , P\si) ■ As in appendix [XI where we follow the 
approach of paper I, it is convenient to first isolate the 
equations of motion in the orbital plane for the variables 
(r, $) and their canonical momenta (P^, P<i)- 
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r = APr 4 




Pr 


= A^l+CPr + D + W^^''^;^ 


$ = A— 


+ ^^"3'^ *(J-L), 




= 



where A,B,C,D arc functions of (r, P^) to be defined 
momentarily. The momentum P$ = L conjugate to $ is 
conserved, by the justification in paper I which survives 
despite the addition of a second spin. In any other basis, 
although L is constant, it is not a momentum conjugate 
to any coordinate. Instead, L should be interpreted in 
terms of the linearly independent coordinates and mo- 
menta appropriate for that basis. The added beauty of 
this non-orthogonal approach is that L is a, canonical mo- 
mentum, namely P$ , not so in the usual spherical coordi- 
nate basis where L is neither a coordinate nor a momen- 
tum and the Hamiltonian angular equations of motion 
are less transparent. 

The four Eqs. (|12p describe motion within the orbital 
plane. The orbital plane itself precesses around the con- 
stant J with variable rate derived in appendix I A 21 to 
be 

\ JxLr^ ) \ J 

(13) 

and P* = L • J = L2 is not a constant. 

Again, using the manipulations of appendixIXl particu- 
larly Eq. (|A27P and the identities J-L = cos 6*^ = P* /P$, 



|J X L| = siuOl, we can write the final term in the $ 
equation of (|12p as 

_ / (Seff-j)-(Seff-L)PWP^ \ 

(14) 

Writing it in this form exploits the dependences on 
the coordinates, conjugate momenta, and the constant 
Seff • L. The one term that clearly remains dependent on 
angles is the term Ses • J. Therefore, when both black 
holes spin, the angular equations will depend on the an- 
gular precession of the orbital plane. 

We saw in paper I a dramatic simplification in the case 
of one effective spin. As follows from the earlier discus- 
sion of the constants of motion, Ses • J would be constant 
if either of the spins vanished and this would remove the 
angular depedence in the above equations. A pair of 
spinning black holes of equal mass is also reducible to a 
system with effectively one spin (see 3p . We continue 
to consider the general case of two misaligned spins for 
arbitrary mass ratios. 

For completeness, and as a complement to paper I, we 
explicitly write out the functions A, B, C, D that were set 
up in as derivatives on the Hamiltonian: 



A=l+l{3v-l)p^-{3 + v)- + 
2 r 



2 (5 - 2O77 - 3rj^) - 2r]' (n • p)^ - + (5 + 8??) ^ 



1 



1 



^ (-5 + 35r; - 70r;2 + 35773) (p^)^ 
1 



3 (-7 + 4277 - 53772 - 577^) (p2)2 + 2 (2 - 3r/) 77^ (n • p)^ p^ + 3 (1 - 77) 77^ (n • p)^ 



27 + 13677 + 10977^) p^ + - (17 + 30?7) 77 (n • p)' 



25 



1 2 335 



64 



23 , 



(15) 



6 



B 



T]{n-p)- + 

r 

-Arf (n • p) - I2rj^ (n • p)' 



1 



3?7 (fi-p)— + 



1 

16 



2 (2 - 37?) if (n • p) (p2)' + 12 (1 - r;) 77^ (fi • p)^ p^ - 30??^ (n • p)^ 



i (17 + 3O77) ?7 (n • p) p2 + 1 (5 + 43r;) 77 (ft • p)' 



r + 



85 

16 " 64' 



- 7?? ) (n ■ p) 



1 



r 

D = -(ft.p)C-l-i ((3 + 77) p2 + 77(ft.pf) 1 + 1- 



(16) 
(17) 



(5 - 20r; - 377') (p2)^ - 2rf (ft • p)^ p^ - 377' (ft • p)^ 



(5 + 877) (p2) + 377 (ft • p)^ 



1 

16 
2 



(-7 + 4277 - 5377^ - 577^) (p)^ + (2 - 377) 77^ (ft • p)' (p^)^ + 3 (1 - 77) 77^ (fi • p)'' p2- 577^ (ft • p 



1 

16 



-27 + 13677 -t- 1097/2) (p2 



2 „2 



16 



(17-f- 307/)?7(ft-p) p 



1 (5 + 4377) 77 (ft -p)^ 



1 



85 3 
16 ^ 64' 



7 



77 77 (fi • p)' 



25 


(1 




V64 




"1 






^4 


8 + 



335 

" Is" 

109 21 

U ^ 32 



23 



TT^ I 77 



(18) 



r 



where fi • p = and p^ = + I? jr^ . Notice that 
A, B,C', D, which come from the non-spinning part of the 
Hamiltonian 0], depend only on {r,Pj.) and constants. 

Useful results can be drawn from a simple observation. 
The radial equations in \12\ have no angular dependence. 
The energy, angular momentum, and radius of spherical 
orbits can be derived from the radial equations alone. 
Therefore we can find spherical orbits simply despite the 
precession of the orbital plane. 

The fact that the two equations in (r, Pr) form a self- 
contained system is a restatement of the fact that the 
Hamiltonian itself can be viewed in a one-dimensional 
effective approach as a function of (r, Pr) and constants. 
It is important to be cautious however when investigat- 
ing the angular motion. The Hamiltonian depends only 
on (r, Pr) and constants in time, yet those constants in 
time have to be carefully varied as functions of the angu- 
lar coordinates and their conjugate momenta in a given 
basis to correctly derive the remaining equations of mo- 
tion. This accounts for the labor in appendix [X] needed 
to derive the ($, P$) and (\1/, P*) equations of motion. 

Still, the simple dependences of the Hamiltonian allows 
us to analyze the spherical orbits as one-dimensional ra- 
dial motion in a simple effective potential. The location 
of the spherical orbits was implicit in paper I to frame 



the distribution of all other orbits. For completeness we 
determine the range of spherical orbits with an eye on 
that companion work. 



II. SPHERICAL ORBITS 

A. Effective Potential for Spinning Black Holes 

Ideally, in an effective potential formulation, the radial 
equation could be cast in the form: 



1 



effective potential = constant 



(19) 



where the effective potential depends only on r and con- 
stants of the motion. Now, the Hamiltonian of Eqs. 
(|4l) does not admit a simple effective potential formula- 
tion since it is a complicated function of p^. We have 
already argued that -ff (r, p, Sgff ) can be written as an ef- 
fective function of {r,Pr) and constants, yet it remains a 
polynomial function of Pr- However, if we only consider 



Kff = H(Pr = 0) 



(20) 



then we have a good representation of a pseudo effective- 
potential at the turning points. We cannot misuse the 
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FIG. 3: An effective potential for two spinning black holes 
ai = a2 = 3/4 of mass ratio -012/ mi = 1/4 for different values 
of the angular momentum. Notice the change in scale between 
panels. Upper; The appearance of the ibso is marked by the 
effective potential touching the line H — 0. Next: As the 
angular momentum decreases, the potential will have both 
stable and unstable spherical orbits. Next: As the angular 
momentum is further decreased there occurs a critical value 
at which the unstable and stable spherical orbits merge at a 
saddle point, the isso. Lower: The last panel shows a differ- 
ence from the Schwarzschild or Kerr stories. At angular mo- 
menta and radii below the occurence of the isso, there occur 
new sets of stable and unstable spherical orbits. These occur 
at radii far below which the approximation can be trusted, 
yet we point out their presence for completeness. 



VcS by trying to interpret motion away from the turning 
points, but it gives a perfectly valid description of the 
behavior at aphelia and periastra as well as on spheri- 
cal orbits. Hereafter we'll shorthand the term "pseudo 
effective-potential" by "effective potential" . 

From the Hamiltonian of Eqs. ([4]), the effective poten- 
tial. 



Vosir,L,Seff ■ £,77) 



(21) 



is a function of orbital parameters (r, L, Sgff • L) and the 
mass ratio. Again, since L and Seff ■ L are constants of 
the motion, for a given (L, Seff • L) and a given mass 
ratio, the potential is a function of r only. 

Fig. [3] shows several snapshots taken of the effective 
potential for a pair of spinning black holes as the mag- 
nitude of L decreases for a given Seff • L value. (For a 
detailed exposition on interpreting effective potentials for 
black hole orbits see Refs. Hi and [l3,[il.) The spher- 
ical orbits are simply the cxtrcma of the potential.^ An 
example of such an orbit was shown in Fig. [T] Although 
this orbit is not generally periodic, it does close in the 
orbital plane as shown in the lower panel of Fig. [TJ 

The top panel of Fig. [3] marks the value of L for which a 
margnially bound, unstable spherical orbit appears. An 
orbit is marginaly bound if its energy H = and it is 
spherical and unstable if it is a maximum of the effective 
potential. The conditions are summarized as 

Kfr(Pr = 0) = 
dV,s 



dr 



= 



< 



(ibso). 



(22) 



although the first two arc sufficient. We call the marg- 
nially bound unstable spherical orbit "ibso" in analogy 
with the innermost unstable circular orbit (ibco) of equa- 
torial orbits. 

For angular momenta below Libso there will be both a 
stable and unstable, energetically bound spherical orbit, 
as in the second snapshot of Fig. [31 until the angular 
momentum gets so low that we reach the third snapshot 
from the top. Here, the unstable and stable spherical 
orbits have merged in a saddle point, coined an innermost 
stable spherical orbit (isso): 



dV, 



cff 



dr 

Qj.2 







(23) 



Orbits with the same angular momentum as a stable spherical or- 
bit but different energy will oscillate between two turning points, 
both of which can be read off the effective potential diagram. 
Again, due to spin precession, for misaligned spins these eccen- 
tric orbits lift out of a plane. Their spectra was shown in paper 
I for a spin/spinless black hole pair. 
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The story plotted out by panels 1-3 of Fig. [3] for 
Lisso < L < iibso qualitatively follows the fully relativis- 
tic Schwarzaschild and Kerr stories as expected ■ 
However, something peculiar then happens in the PN ap- 
proximation at very low values of the angular momentum. 
New stable and unstable spherical orbits can appear as 
shown in the bottom panel of Fig. O Or, for some other 
ranges of parameters, the ibso disappears or the isso dis- 
appears or both disappear. Sometimes these problems 
occur, as in the figure, for radii far below the confidence 
of the PN approximation. We point out these trouble- 
some features in the spirit of full disclosure. More than 



this, the details provide a quantitative testing ground for 




5 


10 


15 


20 


the approximation. 










Despite these oddities at low r where the PN- 
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B. Orbital Parameters for Spherical Orbits 
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For a given black hole pair, that is, a given mass ra- 
tio and Seff • L, all orbits are uniquely specified by their 
[E, L) . Using the effective potential, we can easily gener- 
ate the E and L for spherical orbits and thereby generate 
initial data for them. Initial conditions for spherical or- 
bits were also found in [l^] ■ Damour noticed that 
when only spin-orbit terms are included that the Hamil- 
tonian could be expressed as a radial function. One could 
arrive at this conclusion, as we have in the previous sec- 
tion.^ The constant radius orbits occur at the extrema 
of H(r,Pr), in the same spirit as an effective potential 
method. 

From the vantage point of the effective potential, 
spherical orbits satisfy the condition 



FIG. 4: {m2/m\ = 10 ^,5*011 = 0). Upper: Angular momen- 
tum vs rs- Lower: Energy vs r^. 



dr 



= 



(24) 



treating L and Seff • L as constants. We could also take 
the vantage point of the equations of motion. Since, P,. = 
forces B = 0, the condition f = AP,- + 5 = can be 
thought of as synonymous with the condition that Pr = 
0. The constant radius condition is thus the requirement 
that 







(25) 



in Eqs. (HH) and Eq. ^ is equivalent to Eq. ([M)) . 

The roots of the 8th-ordcr equation in L, Eq. ((25)) . give 
the angular momenta of the spherical orbits as a function 
of spherical radius, Lg[rg), where we use a subscript s to 
denote a quantity evaluated at a spherical orbit. (When 
there is no spin, the condition reduces to a quartic in 
with only two of the four roots real.) Piecing together 
the real roots we find L^'s such as the one in Fig. U) 
Although the upper branch grows very quickly in Lg, 
these values rapidly become physically unreachable since 
it would require angular velocities greater than the speed 
of light to be that high up on the upper branch. One can 
think of [Ls/ts) < 1 as a crude marker of physically 
allowed values. To find the energy of spherical orbits, 
Es{rs), we simply plug Ls{rs) into the Hamiltonian when 
P,- = 0. The energy plot is also shown in Fig. 01 

There are several things to notice about the and 
Eg plots, for which wc chose values to illustrate the PN 
approximation to Schwarzschild (m2/mi ~ 10~^,S'cff = 
0).^ The large values of rg correspond to stable spherical 
orbits. (Since spin is zero, these constant radius orbits 
are actually circular equatorial but we'll keep the lan- 



However, in another basis such as the usual basis for spherical 
coordinates used in [loll , projection of the vector equations of 
motion will give equations of motion in component from that 
continue to depend on angles even when only one black hole 
spins, unlike the orbital basis of i|A 31 



Since the radial equation depends only on the combination Seff • 
L, Fig.|4]should be equally valid for a non-zero effective spin that 
is orthogonal L. 
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guagc more general.) When hits a minimum, we have 
found the isso - for no L < Ljsso ^ire there spherical orbits. 
To the left of that minimum arc the unstable spherical 
orbits. A true peculiarity of the figure is the fact that the 
radii of the unstable spherical orbits begin to move out to 
larger r. This is simply a flaw in the PN approximation 
and does not occur in the Schwarzschild system. In the 
fully relativistic system the unstable spherical orbits al- 
ways move to smaller radii than the isco, hence the ibco 
is really innermost, earning its name. Here, the ibso is 
not actually innermost - due to the poor quality of the 
approximation - although it remains the highest energy 
bound spherical orbit when it exists. The ibso cannot be 
read off of Fig. [4] although it can be found simply as the 
coincident of the roots of Eqs. ((22|) . 

Figs. [4] are for a non-spinning extreme-mass-ratio bi- 
nary and are therefore valid as an approximation to 
Schwarzschild. The details of these figures will be use- 
ful for a future test of the PN expansion. Here we note 
that Libso ~ 4.69, which is about 17% higher than the 
Schwarzschild value of 4 while iisso ~ 3.75, which is 
about 8% higher than the Schwarzschild value of vT2. 
The energy of the ibso is designed to be zero so is not 
informative but the energy of the isso is i?isso ~ —0.0452, 
which is about 21% less negative, that is less energeti- 
cally bound, than the Schwarzschild case (2-\/2/3) — 1. 
Due to the approximate nature of the expansion it is not 
necessary to take these comparisons to heart, but they 
indicate how the spherical orbits and the periodic spectra 
could facillitate a test of the PN expansion. For a com- 
parison of the isso in different PN approaches including 
the resummed Kerr-like effective-one-body approach see 

MM. 
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FIG. 5: All black hole pairs represented have Seff ■ L = 
0.35355. Upper: Angular momentum vs r for the ibso and 
isso for different mass ratios. The upper point is always the 
ibso for a given symbol while the lower point with the same 
symbol is always the isso. The key lists the different (m2/mi). 
Lower; Energy vs r. 



C. Dependence of Binding Energies on mass ratios 
and spin 

For completeness, we can see how the ibso and the isso 
vary as the mass ratio and spins of the black holes are 
varied; that is, as their mass ratio and spins are varied. 
For one, since the ibso and isso frame the distribution of 
orbits, they define the ranges of E and L values for all 
other orbits in the strong-field. For another, the energy 
at the isso gives an estimate of the energy emitted on 
quasi-circular inspiral up to the transition to plunge. A 
larger binding energy at the isso could also mean a larger 
signal at final coalescence so these variations attest to 
various levels of detcctability. We will discuss the binding 
energy of the isso in this section. In the next section we 
will consider the transition to plunge for eccentric orbits. 

A black hole pair is specified by its mass ratio, 7712 /mi, 
and its spins through the particular combination Seff • L. 
The Hamiltonian, and the radial equations, depend only 
on these two combinations. We will therefore consider 
the variations in the isso and ibso for black hole pairs 
distinguished only by their (m2/mi, Seff • L) values. It is 
important to realize that there is a great deal of degener- 



acy among pairs. The ibso and isso values (their energy, 
angular momenta, and radial values) are identical for two 
physically distinct black hole pairs. For instance, a black 
hole with mass ratio m2/mi = 1/3 and Seff -L = 0.35355 
could be a black hole with initial values ai = 1/4, a2 = 0, 
and the spin of the heavier black hole aligned with the 
initial orbital angular momentum. However, this is not 
the only combination of spin amplitudes and angles that 
will give the combination Seff • L — 0.35355. While dif- 
ferent black hole pairs can give degenerate isso and ibso 
values, they will be physically distinguishable through 
their angular motion. 

In Fig. [5] the L of the ibso and isso is plotted in the up- 
per panel and the E of the isso and of the ibso are plotted 
in the lower panel. Qualitative conclusions can be drawn 
from these figures. Wc notice that as the mass ratio is 
increased towards 1 , the radius of both the isso and ibso 
decrease, although the isso moves in faster. Therefore 
the isso is pushed to larger binding energies as the mass 
ratio is increased towards 1. Because the Hamiltonian 
is a high-order polynomial in r, there can be more than 
one marginally bound orbit and more than one saddle 
point for a given (m2/mi, Seff • L) pair, as demonstrated 
in the lowest panel of Fig. [3] The second occurence of a 
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FIG. 6: Upper: Angular momentum vs r for the ibso and isso 
for fixed mass ratio mi /mi = 1/3 but varying S^ff ■ L. The 
upper point is always the ibso for a given symbol while the 
lower point with the same symbol is always the isso. The key 
lists Seff • L. Lower: Energy vs r. 



marginally bound orbit and/or saddle point appears in 
the vicinity of r 1 where the approximation is unin- 
terpretable. (There may even be third occurences.) It 
is unclear if there is any physical content to these other 
stable and unstable spherical orbits. Fig.[5]plots only the 
ibso/isso pair for r values > 2. 

For the value of Seff • L ~ 0.35355 used in the figure, 
either the ibso or the isso disappears (or both disappear) 
as m2 approaches mi. There may still be very small 
radii (r ^ 1) ibso's and/or isso's, but the sensible ones 
disappear. This peculiarity is probably an artefact of the 
approximation, a point we return to momentarily. 

Fig.[H]fixes the mass ratio at milm\ = 1/3 and varies 
Seff • L. Increasing Seff • L has the same effect of pushing 
the isso to smaller separations and therefore to larger 
binding energies, although again the isso moves in faster 
- discounting any marginally bound spherical orbit or 
saddle points that occur in the vicinity of r ^ 1. (In 
fact, as Fig. [5] shows, at some point the isso actually 
occurs at a smaller radius than the ibso.) So, all other 
factors being equal, spins anti-aligned with the orbital 
angular momentum push the isso out to larger radii and 
smaller binding energies while aligned spins pull the isso 
into smaller radii and larger binding energies. For the 



mass ratio of this figure, the ibso or the isso actually 
vanishes (or both vanish) as Seff • L is increased much 
beyond the values shown. 

These trends are consistent with those for spherical or- 
bits discussed in Ref. [l^l- For the equal mass case with 
spins aligned or anti-aligned with the orbital angular mo- 
mentum, the authors of that reference also remarked on 
the absense of an isso (called a last stable spherical or- 
bit (Isso) in their lexicon). Indeed, they used this fail- 
ing to argue that the PN expansion could not be used 
to study the transition from inspirl to plunge and advo- 
cated instead the seemingly more reliable effective-one- 
body (FOB) approach [i,!!!,!!^]. It would be interesting 
to extend to the EOB method the investigation of the 
zoom-whirl orbits of paper I Q and the homoclinic limit 
of the zoom-whirls that we turn to in the next section. 
We leave that to a future work and continue to use the 
3PN Hamiltonian as an example of our general method. 

The disappearance of the ibso accompanies the disap- 
pearance of all unstable spherical orbits. Once this hap- 
pens, there can be no isso since the isso is really the point 
of merger of the unstable and stable spherical orbits.^ 
As is already known, in the absence of spin there are no 
bound unstable circular orbits at 2PN (23l. [29l. [3l]. [3^. [33| . 
At 3PN there are no bound unstable circular orbits for 
mass ratios bigger than about m2/mi ~ 1/2. (See also 
pn. Is^l-) The absence of a bound unstable circular or- 
bit is clearly a shortcoming of the approximation since we 
know that the Schwarzschild spacetime possess an unsta- 
ble circular orbit as a refiection of its high non-linearity. 
Furthermore, the unstable circular orbits are present in 
fully rclativistic treatments as the equal mass numerical 
investigation of Ref. [2^ shows. Therefore, the ibso and 
isso should emerge for m2 mi at higher orders. In- 
cidentally, their disappearance at 3PN-order implies the 
expansion is very likely approximating the dynamics as 
more stable than it really is and therefore less vulnera- 
ble to chaos than it really is for these comparable mass 
binaries. 

We have already warned caution to take the trends as 
qualitative indicators and not to invest too much in the 
numbers due to pressures on the PN expansion at such 
large values of (mi + m2)/r. Afterall, the PN expansion 
is an expansion in small (mi + m2)/r and will naturally 
begin to faulter for small r. 

We have focused on the binding energy of the isso pri- 
marily to fit into the wider conversation that has focused 
on quasi-circular inspiral. However, the eccentric binaries 
formed by tidal capture in dense regions will not transi- 
tion from inspiral to plunge through the isso. Rather 
they will transition through the eccentric separatrix be- 
tween bound and plunging orbits. We investigate that 
separatrix briefly in the final section. 



This is not the only reason the isso disappears. Sometimes the 
isso disappears because the potential simply never flattens out. 
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III. HOMOCLINIC ORBITS - THE -""50 
SEPARATRIX BETWEEN BOUND AND 
PLUNGING ORBITS 

-0.055 



It is worthwhile to mention another important kind 
of orbit that occurs in our dynamical system, the ho- 
moclinic orbit 0, [l2| ■ Homoclinic orbits arc intriguing 
for several reasons, not least of which is that they mark 
the orbits through which the transition from inspiral to 
plunge should occur. In fact, the isso itself, the transi- 
tion point for quasi-circular orbits, is a zero eccentricity 
homoclinic orbit [13, 0, ■ We make the connection be- 
tween the energetically bound, unstable spherical orbits 
and the homoclinic orbit explicit in this final section. 

Formally, homoclinic orbits are defined as trajectories 
that asymptote to the same hyperbolic invariant set in 
the infinite future as in the inifinite past. In these black 
hole settings, the role of the hyperbolic invariant set is 
played by the energetically bound, unstable spherical or- 
bits. Although in the lexicon of black hole physics these 
orbits have been coined "unstable", they are strictly 
speaking hyperoblic, which is to say they possess both 
a stable eigendirection and an unstable eigendirection 
under linear perturbations. And, the eigendirections lie 
along the homoclinic orbit in the local neighborhood of 
the unstable circle they approach. Although we won't 
demonstrate that line up here, the point was emphasized 
in detail in Refs. [H, [iJI for Kerr equatorial dynamics. 

The stability exponents can be found by linearizing in 
small perturbations around E as. (fT2l) . This was done for 
equatorial Kerr orbits in Ref . [1J| . Although we will not 
write out the explicit procedure here, we mention that 
the spherical orbits have radial eigenvalues that come in 
plus/minus pairs, as they must in a Hamiltonian system. 
The radial eigenvalues arc real for the unstable spherical 
orbits and imaginary for the stable spherical orbits. The 
isso occurs at the merger of the eigenvalues at zero. A 
direct computation of the stability exponents around cir- 
cular orbits confirms that the stable spherical orbits and 
the unstable spherical orbits arc distributed around the 
isso as Fig. m shows. 

Through the phase space analysis we have shown that 
the energetically bound, unstable circular orbits arc ac- 
tually hyperbolic - they have a positive stability expo- 
nent as well as a negative stability exponent. We could 
compute the eigenvectors and show they lie along the ho- 
moclinic orbit in the local neighborhood of the unstable 
circle as we did for Kerr in Ref. [l4| . However, for our 
purposes it is sufficient and illuminating to consider a 
physical space picture. 

We can identify the separatrix - i.e. the homoclinic 
orbit - in an effective-potential picture. In particular, 
consider a binary with mass ratio m2/mi = 1/4 and the 
heavier black hole spins with amplitude ai = 1/2 offset 
from L by 7r/4 while the lighter black hole is nonspinning. 
Although, again, any equivalent combination of Sgff ■ L 
is described by this same figure. The unstable spherical 
orbit, r„, at the maximum of Vcs in Fig. [7] is drawn in 
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FIG. 7: An effective-potential for 1x12/ m\ — 1/4, with the 
spin of the heavier black hole displaced from L by 7r/4 and 
amplitude ai = 1/2 while the lighter black hole has no spin. 
The straight line is the energy of the unstable spherical orbit. 
It is also the energy at another, larger turning point Va ~ 10, 
which identifies the apaastron of the homoclinc orbit. 

physical space in Fig. [S] Although the orbit is a closed 
circle in the orbital plane, it fills out a band on a sphere in 
three dimensions. Because of numerical instability near 
this orbit, we only show a few windings. 

The energy of this orbit, Es{ru), is indicated by a 
straight line across the potential of Fig. [T] Note that 
this energy touches the potential at the unstable radius 
r„ and at some larger radius, roughly r ~ 10. This larger 
radius is the apastron of an orbit. If the two black holes 
arc released from rest at an initial separation in center-of- 
mass coordinates equal to this apastron, their orbit will 
roll down the potential (although the shape changes when 
Pr 7^ 0) and then climb back up the other side asymptot- 
ically approaching the spherical orbit at the top of the 
hill. By definition, this is a homoclinic orbit (Fig. (5]). To 
our knowledge it is the first of its kind to be found out 
of the equatorial plane [1, [l^, [3l • 

The orbit winds around the center of mass an infinite 
number of times as it asymptotically approaches the un- 
stable spherical orbit. Although not strictly periodic - 
the homoclinic orbit never returns to apastron - it will be 
significant for the periodic tables [3, [1] as a maximum en- 
ergy orbit for a given L in the strong-field regime (isl . [l^ . 



IV. CONCLUSIONS 

This paper, as the second in a scries, provides the en- 
ergetic frame in which the periodic tables of paper I were 
set 01 . Although in support of paper Fs goals, the anal- 
ysis of spherical orbits could be relevant to additional 
tests of the FN expansion for spinning black hole pairs 
and could have a place in the disucssion of initial values 
for numerical relativity. Additionally, we find the non- 
equatorial homoclinic orbits that whirl an infinite num- 
ber of times as they asymptote to the unstable spherical 
orbit. The homoclinic spearatrix is important as defining 
the transition to plunge for all orbits, including eccentric 
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FIG. 8: The unstable spherical orbit that is the maximum of Kff for Fig. [7] Unlike the effective potential, the details of the 
full orbit do depend on the specific combination Seff ■ L. Left: As viewed in three dimensions. Right: As viewed in the orbital 
plane. 
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FIG. 9: The homoclinic orbit for Fig. [7] approaching the unstable spherical orbit. Left: As viewed in 3d. Right: As viewed in 
the orbital plane. Because of numerical instability near the highly unstable constant radius orbit, we only show a few windings. 



and precessing orbits. It would be interesting to extend 
this study to the EOB method 0, [H, HO] in a future 
work. 

In closing, we comment on an intriguing implication of 
the set of spherical orbits for spinning black hole pairs. 
In this work, we restricted ourselves to spin-orbit cou- 
pling and although we allowed both black holes to spin, 
we found that the spherical orbits constrain the range of 
allowed bound orbits in the following sense. For a given 
angular momentum, spin initial conditions, and mass ra- 
tio, the stable spherical orbit is the lowest energy geodesic 
and the unstable spherical orbit is the highest energy 
orbit in the strong-field - barring the failures of the ap- 
proximation at these close separations.^ As we showed in 
paper I Q , between these two spherical orbits lies an infi- 
nite set of orbits that are closed in the orbital plane. The 
periodic set corresponds to a subset of the rationals, with 



® We actually consider the emergence of an ibco to define the 
strong-field. For the equal mass cases that resist the develop- 
ment of an ibco, it is as if the approximation is not effective 
enough to enter the strong field. 



the rational identifying a given orbit increasing monoton- 
ically between the stable spherical orbit and the unstable 
spherical orbit. The homoclinic orbit is the infinite whirl 
limit of the periodic set and would be the final entry in a 
periodic table of orbits corresponding to the infinite limit 
of the rationals. 



This pattern of a periodic set framed by constant ra- 
dius orbits and limiting to the homoclinic is consistent 
with a pict ure that has emerged for Kerr black hole or- 
bits [1,111, [T^ . The consistency of the picture for two 
spining comparable mass black holes with the Kerr case 
is precisely what is surprising, or at least intriguing. The 
geo desies in a Kerr spacetime are known to be integrable 
[35j . There are enough constants of the motion to restrict 
trajectories to regular tori and prohibit chaotic mixing. 
As Poincare intuited, the structure of the periodic or- 
bits encodes the entire dynamics and the regularity of 
the system is in fact reflected in the regularity of the pe- 
riodic spectrum. The simplicity of the spherical orbits 
and the periodic set they frame suggests that even when 
both black holes spin and are of comparable mass, there 
is no chaos - at least not in physically plausible regimes 
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-if only spin-orbit coupling is included. 

Put another way, homoelinie orbits are also a sign of 
non-linearity. They mark the intersection of the sta- 
ble and unstable manifolds of a hyperbolic invariant set. 
They are the precursor to chaos in the sense that un- 
der perturbation, the homoelinie orbit breaks up into a 
homoelinie tangle and will be the locus of a fractal set 
of orbits [3, [3^. The fractal set is sometimes refered 
to as a strange repellor and is the analog for conserva- 
tive systems of strange attractors in dissipative systems 

Systems with a regular set of periodic orbits that cul- 
minate in a homoelinie limit are not chaotic. However, 
the spinning pairs are vulnerable to chaos as evidenced by 
their very possesion of a homoelinie orbit. Indeed chaos 
has by now been well confirmed in the form of a fractal set 
when spin-spin coupling is included [^, 
As suspected in Ref. [43, our work suggests that the 
emergence of chaos must be directly tracable to the spin- 
spin coupling. We conjecture that the transition to chaos 
could be witnessed through the destruction of the cor- 
resondence of the periodic set with the rationals when 
spin-spin coupling is turned on. The additional proces- 
sional effects of spin-spin coupling, we suggest, must de- 
stroy the homoelinie orbit, replace it with a homoelinie 
tangle - a fractal set of orbits - and induce chaotic scat- 
tering among geodesies in the vicinity. 
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It is possible that for S^ff ■ L much larger than black hole physical 
values chaos could develop. Afterall, one of the constants of 
motion Pip has been lost with the inclusion of spin-orbit coupling 
opening the door for chaos. 
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APPENDIX A: PROJECTION THE EQUATIONS 
OF MOTION ONTO THE NON-ORTHOGONAL 
ORBITAL BASIS 

1. The Orbital Plane Equations 

The four equations of motion in the orbital plane are 
obtained by projecting Hamilton's equations onto the ba- 
sis vectors, as is done in celestial mechanics. For now, 
consider only the projections onto the orbital basis vec- 
tors to generate the four equations, 

. . OH ^ 



r • $ 
p • n 
p • $ 



dH 

dH 
dr 
dH 

dr 



n 



(Al) 



To break down the LHS involves 



Rn + Prh — -^r^ + — $ 



(A2) 



We will need projections of n and $ along n and Now, 
since n • fi = 1, it follows that n • n = and by the same 

reasoning $ • ^ = 0. Also, by orthogonality. 



n . ^ = ^ n 



(A3) 



To obtain the final dot product above we expand the 
basis vectors (n, '4') in terms of an intermediate basis 
(X, Y) that spans the orbital plane, and then expanding 
(X, Y) in the Cartesian basis. We proceed as we did in 
paper I and define the intersection of the orbital plane 
with the equatorial plane: 



J X L J X L 

JxL sin^L 



(A4) 



where cos6'l = L • J. The vector orthogonal to X that 
lies in the orbital plane is 



Y = L X X 



(A5) 



This intermediate orbital basis will be useful in the ma- 
nipulations that follow. In terms of Cartesian compo- 
nents defined with k = J and i, j spanning the equatorial 
plane, we can expand 

X = cos *i + sin 'J'j 

Y = sin6'y(-sin\E'i + cosvI'j) + cos6'yk , (A6) 



where cos Oy = Y • J. Since Y is always orthogonal to L, 
again by construction, this is not really a new angle but 
can be recast a.s Oy = 'n'/'i — 9^. 

Our non-orthogonal basis can then be expanded as 



Using 



X 

4f 



n 



-'i'X 



cos $X + sin $Y 

- sin $X + cos <i>Y 

— sin 'fi + cos vj/j 



(A7) 

(A8) 
(A9) 



Y = - sin 6iy*X + cos 6'y6'y*- sin 6iy6lyk (AlO) 

From all of the above relations we obtain for use in the 
projections 

A-n = (All) 

* • * = (A12) 

= ^ + 4/sm6y = ^ + i'cos6L (A13) 

* • n = -n • $ . (A14) 

Conveniently, these are the same projections we found in 
paper I for the case in which only one black hole spins 
and 6y =eL = 0. 

Now we can derive the equations of motion in the 
(r, ^P) coordinates. We use the equations we con- 
structed in paper I 0] 



r = yip + Bn 



Seff X r 



p = Cp + Dn H h 3 7 — n 



(A15) 



where A,B,C,D are given by Eqs. p8|) . With the pro- 
jections (Eqs. (jAip ). (|A2p . and the above vector relations 
we have the radial equation from r • n in (jAlSp : 



APr + B 



The $ equation follows from 

dH - 
r . $ = __ $ 

dp 



\ / r r 



Look at 

(Seff xr).* = r(Seff -L 

The $ equation is then 



• L ■ „ Seff • L 



(A16) 

(A17) 
.(A18) 

(A19) 
(A20) 



where S • L is constant. 
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The two conjugate momenta equations are next. We 
start with Pr'. 



p-n = P,, + * cos 6*^ 



(A21) 



'eff 



eff 



= CPr+D + 2 

where we have used that 

(p X Seff) • n = 

r 

Notice if we use Eq. (jA20[) . we have 

P, = AK+CPr + D + 3^^^ 



(A22) 



and last 



p • l> = Pr{^ + ^1 cos6l) 



(A23) 



where we have used that 

(p X Seff) • * = Seff • (* X p) = -P,Seff • L (A24) 

Notice if we use Eq. (|A20|) . we have a cancellation and 

which confirms a true statement but docs not provide any 
new equation of motion. The final equation of motion is 
simply Pcj) = 0. All four equations in the orbital basis 
are compiled in the boxed Eqs. ([l3 



2. The Precession of the Plane 

The plane precesses in the direction ^ at a rate 5*, 
which can be computed from the first of Eqs. (|A10[) : 

i = 

We can isolate by projecting along ^f, 

X • * = * . (A25) 

We take the time derivative of Eq. (jA4p and use the con- 
stancy of J and the precession equation for L from Eq. 
(H) to find 



* = 



j X I Seff X L 



J X L r-' 



(A26) 



Notice that the term that would have been proportional 
to 0i is killed since it is also proportional to X • ^ = 0. 
With some vector manipulations, including the general 
rule Ax (Bx C) = B(A • C) - C(A • B), applied to both 
the term in parantheses and to ^ = J x X with X given 
by Eq. (|A4|) . this can be reduced to 



Seff • (j - L(J • L) 

sin 6\r^ 



(A27) 



Going in the other direction with the general rule, B(A • 
C) — C(A-B) = Ax (Bx C), we can write the right-hand- 
side as a triple cross product and identify the particularly 
compact form 



5 eff 



sin O^r^ 



As in paper I, 



= = L • J 



(A28) 



(A29) 



Unlike paper I, P* , which can also be expressed as P^/ = 
L cos , is not conserved when both black holes spin and 
precess. 



3. One Effective Spin 

The equations of motion simplify considerably if there 
is only one effective spin, such as the case of only black 
hole spinning 



f = APr 










Pr 


^A^l ^P,. + D + 36,^'-^^ 












- Ol^ , 


P$ 


= 



The orbital plane precesses with frequency 

J 



p* 



0. 



(A30) 



Consequently, the equations of motion above are inde- 
pendent of angles. In paper I, we used these purely radial 
equations to study several features of the dynamical sys- 
tem, such as a periodic table that defined the spectrum 
of black hole orbits. 

The same simplification can be effected when the black 
holes are of equal mass mi = m2- Then what we really 
mean by Si is Si ^ Si + S2. 



16 



D 



S. F. Portegeis-Zwart and a.-p. S L McMillan. 
L. Wen, Ap. J. 598, 419 (2003). 
R. M. O'Leary, B. Kocsis, and A. Loeb, (2008). 
J. Levin, R. O'Reilly, and E. Copeland, Phys. Rev 
62, 024023 (2000). 

J. Levin, Phys. Rev. Lett. 84, 3515 (2000). 

J. Levin, Phys. Rev. D 67, 044013 (2003). 

J. Levin and R. Grossman, gr-qc/08093838 (2008). 

J. Levin and G. Perez-Giz, Phys. Rev. D 77, 103005 

(2008). 

T. Damour, Phys. Rev. D 64, 124013 (2001). 

A. Buonanno, Y. Chen, and T. Damour, Phys. Rev. D 

74, 104005 (2006). 

U. Sperhake et al., Phys. Rev. D78, 064069 (2008). 
Bombelli and Calzetta, Class and Quant. Grav. 9, 2573 
(1992). 

J. Levin and G. Perez-Giz, Homoclinic Orbits around 
Spinning Black Holes I: Exact Solution for the Kerr Sep- 
aratrix, 2008. 

G. Perez-Giz and J. Levin, Homoclinic Orbits around 
Spinning Black Holes II: The Phase Space Portrait, 2008. 

F. Pretorius and D. Khurana, Class. Quant. Grav. 24, 
S83 (2007). 

G. Schaefer, Annals Phys. 161 
T. Damour and N. Deruelle, C. 
537 (1981). 

T. Damour and G. Schaefer 
(1988). 

P. Jaranowski and G. Schaefer, 
029902 (2001). 

T. Damour, P. Jaranowski, and G 
B 513, 147 (2001). 
T. Damour, P. Jaranowski, and G. 
D 62, 084011 (2000). 

A. Gopakumar and C. Konigsdorffer, Phys. Rev. D 72 



, 81 (1985). 
R. Acad. Sci. 



Paris 293, 



Nuovo Cim. BlOl, 127 
erratum-ibib.d 63, 
. Schafer, Phys. Lett. 
Schaefer, Phys. Rev. 



[23' 
[24' 
[25^ 
[26; 

[27 
[28' 
[29 

[30 



[31 

[32' 
[33; 

[34' 
[35' 
[36 

[37 

[38; 

[39 

[40 
[41 
[42' 
[43 



121501 (2005). 

Schafer and Wex, Phys. Lett. A 174 (1993). 

N. Wex and S. Kopeikin, astro-ph/9811052 (1998). 

B. Gong, astro-ph/0401152 (2004). 

C. Konigsdoerffer and A. Gopakumar, Phys. Rev. D 71, 
024039 (2005). 

R. Wald, General Relativity, 1984. 

F. Pretorius, Class. Quant. Grav. 23 (2006). 

A. Buonanno and T. Damour, Phys. Rev. D 59, 084006 

(1999). 

A. Buonanno et al.. Toward faithful templates for non- 
spinning binary black holes using the effective-one-body 
approach, 2007. 

L. E. Kidder, C. M. Will, and A. G. Wiseman, Phys. 
Rev. D 47, 4183 (1993). 

N. Wex and G. Schafer, Class, and Quantum Grav. . 
T. Damour, B. R. Iyer, and B. S. Sathyaprakash, Phys. 
Rev. D57, 885 (1998). 

L. Blanchet, Phys. Rev. D65, 124009 (2002). 

B. Carter, Phys. Rev. 174, 1559 (1968). 

E. Ott, Chaos in Dynamical Systems, Cambridge Uni- 
versity Press, 2002. 

N. J. Cornish, C. P. Dettmann, and N. E. Frankel, Phys. 
Rev. D 50, 618 (1994). 

C. Dettmann, N. Frankel, and N. Cornish, Phys. Rev. 

D. 50, 618 (1994). 

N. J. Cornish and J. J. Levin, Phys. Rev. Lett. 78, 998 
(1997). 

J. Levin, Phys. Rev. D 74, 124027 (2006). 

X. Wu and Y. Xie, Phys. Rev. D 76, 124004 (2007). 

X. Wu and Y. Xie, Phys. Rev. D 77, 103012 (2008). 

M. D. Hartl and A. Buonanno, Phys. Rev. D 71, 024027 

(2005). 



